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THE  ARMY  NEW  PERSONNEL  SYSTEM  EVALUATION  MODEL 


Introduction 

As  a  result  of  the  changes  in  the  funding  levels  of  the 
Department  of  Defense,  the  Army  personnel  system  is  currently 
undergoing  significant  change.  The  overseas  manpower  requirement 
will  be  reduced  and  the  rotation  policies  will  need  to  be 
evaluated.  The  permanent  change  of  station  (PCS)  policies  need  to 
be  reexamined. 

This  report  gives  a  survey  of  techniques  that  have  been 
used  in  personnel  policy  analysis  along  with  a  prototype 
mathematical  model.  The  prototype  model  demonstrates  the  key 
features  needed  in  a  production  model,  along  with  an  illustration 
of  how  such  a  production  model  could  be  developed.  We  also  give 
the  model  generator  code. 


The  Survey 

Manpower  planning  is  concerned  with  the  allocation  of  the 
right  number  of  personnel  to  different  tasks  in  order  to  achieve 
short  and  long  term  goals  of  an  organization  without  violating 
organizational  policies.  The  extensive  use  of  the  term  manpower 
planning  started  after  World  War  II  when  the  U.S.  Navy  began 
reorganization  of  its  technical  and  managerial  manpower.  Manpower 
planning  models  have  been  used  for  government  and  public 
agencies  and  extensively  in  the  military  due  to  the  importance 
of  national  defense  planning  and  budget  issues.  Accurate 
forecasts  and  eminent  knowledge  of  market  trends  are  necessary 
for  a  successful  development  of  a  manpower  system. 

Markovian  models  have  been  used  to  estimate  grade-wise 
distribution  of  future  manpower  as  well  as  to  maintain  new 
recruitment  and  promotion  or  firing  policies  (see  Vajda  [1976], 
Abodunde  and  McClean  [1980],  Zanakis  and  Maret  [1980], 

Bartholomew  [1982],  Edwards  [1983],  and  Raghavendra  [1991]).  A 
combination  of  Markov  models  and  linear  programming  is  used  to 
obtain  optimal  policies  taking  into  account  not  only  the  manpower 
requirements  but  also  costs  and  conflicting  objectives  (see  Young 
and  Abodunde  [1979],  and  Zanakis  and  Maret  [1981]).  Dynamic 
programming  techniques  and  optimal  control  theory  may  also  be 
used  in  manpower  planning  but  are  not  commonly  used  (see  Edwards 
[1983] )  . 

Vajda  [1976]  and  Edwards  [1983]  model  personnel  movement 
within  an  organization  as  a  Markov  chain.  Vassiliou  [1976] 
develops  a  Markov  model  to  analyze  wastage,  departures,  from  an 
organization.  Young  and  Abodunde  [1979]  represent  internal  and 
external  transitions  of  personnel  within  a  graded  organization  as 
a  discrete  time  deterministic  model  and  present  a  linear 
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programming  formulation  of  the  problem  to  obtain  optimal  long 
term  recruitment  policy.  Abodunde  and  McClean  [1980]  analyze  the 
effects  of  possion  recruitment  on  a  stochastic  version  of  the 
Young  and  Abodunde  model.  Zanakis  and  Maret  [1981]  discuss 
different  stages  of  development  and  validation  of  a  Markov  model, 
and  combine  a  Markov  model  and  linear  goal  programming  for 
manpower  planning.  Edwards  [1983]  conducts  a  survey  of  manpower 
planning  applications  in  the  U.K.  and  other  European  countries. 
The  survey  includes  different  software  available  for  manpower 
planning,  the  models  and  the  techniques  used  in  their  development 
as  well  as  requirements  for  using  this  software.  Hall  [1986] 
proposed  an  aggregate  manpower  planning  system  using  graphical 
models.  Smith  and  Bartholomew  [1988]  present  a  historical  review 
that  traces  the  growth  of  manpower  planning  in  the  U.K.  since 
World  War  II.  Raghavendra  [1991]  uses  a  Markov  chain  model  to 
estimate  the  probability  matrix  for  promotions  in  an 
organization.  He  then  derives  the  formulas  for  obtaining  a 
minimum  level  of  seniority  and  performance  rate  required  for 
promotion. 

The  needs  for  manipulating  a  multi-characteristic  pool  of 
personnel  over  a  long  planning  horizon  has  led  to  a  number  of 
linear  programming  based  approaches  for  evaluating  the 
alternative  solutions.  The  linear  goal  programming  model  (LGPM) 
provided  an  indispensable  aid  for  managing  the  flow  of  personnel 
in  a  manner  that  best  meets  the  desired  manning  structure  and 
quality  over  the  planning  years. 

Several  LGPM  based  pragmatic  systems  have  been  developed  and 
extensively  exploited  in  the  US  Army  and  Navy.  Two  such  systems 
are  the  Accession  Supply  Costing  and  Requirement  (ASCAR)  and 
Enlisted  Loss  Inventory  Model-Computation  of  Manpower  Programs 
Using  Linear  Programming  (ELIM-COMPLIP)  models  (see  Collins,  Gass 
and  Rosendahl  [1983],  and  Holz  and  Worth  [1980]).  The  primary  use 
of  these  systems  is  to  evaluate  the  impact  of  changing  the 
manpower  policy,  such  as  variations  in  the  promotion  and 
separation  rates  and  desired  force  compositions,  on  the  overall 
performance  and  man-strength.  Also,  they  have  been  used  to 
provide  the  capability  to  determine  the  number  of  new  recruits 
with  certain  qualifications  to  meet  future  needs,  to  improve  the 
ability  to  forecast  the  budgets  required  to  cover  certain  future 
manning  levels,  and  to  establish  the  ability  to  significantly 
increase  or  reduce  the  man-strength  in  a  short  period  of  time  to 
meet  critical  situations  such  as  Vietnam  and  the  Gulf  war. 

Goal  priorities  are  very  decisive  in  affecting  the  solution 
profile.  These  priorities  can  be  incorporated  in  the  model  by 
associating  suitable  weights  with  the  goal  deviation  variables. 
Determining  these  weights  directly  is  not  feasible,  Gass  [1986] 
developed  an  automated  system  for  calculating  composite  weights, 
that  properly  reflect  priorities  of  composite  goals,  using 
Saaty's  analytical  hierarchy  process  (AHP) . 
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In  the  following  section  we  present  a  summary  of  some  of  the 
Markovian  approaches  to  the  manpower  planning  problem,  followed 
by  a  summary  of  linear  goal  programming  approaches.  A  combined 
approach  is  presented  in  the  last  section. 

Markov  Approaches 

Markov  Approaches 

Manpower  planning  consists  of  demand  prediction,  supply 
prediction,  and  development  of  policies  to  reconcile  any 
difference  between  supply  and  demand.  Future  demands  can  be 
predicted  by  simple  extrapolation,  regression,  factor  analysis, 
work  study,  or  other  forecasting  techniques.  To  predict  the 
future  supply  one  may  represent  the  manpower  system  as  stocks  and 
flows,  where  stocks  are  number  of  employees  in  each  category  and 
flows  are  movement  of  personnel  into,  within,  or  out  of  the 
system.  While  most  flows  are  under  complete  managerial  control, 
voluntary  leaves  cannot  be  totally  controlled  by  management.  An 
indicator  of  the  voluntary  leave  survival  function  may  be 
obtained  using  cohort  or  census  analysis.  If  the  future  demand 
has  been  predicted  and  stocks  and  flows  determined,  possible 
shortage  or  excess  can  easily  be  obtained  (see  Edwards,  1983) . 

A  mathematical  statement  of  the  problem  is  made  possible  by 
classifying  the  personnel  into  different  groups  by  their  sex, 
race,  level  of  employment,  education,  etc.,  and  then  representing 
these  groups  as  different  states  of  a  Markov  model.  The 
transition  probability  matrix  (TPM)  represents  the  movements 
within  the  organization  such  as  promotions  and  transfers  or  any 
rnovem®  it  of  the  personnel  from,  to,  or  within  an  organization. 
States  such  as  retirement,  voluntary  and  involuntary  wastage  are 
absorbing  states,  while  others  are  transient  states.  To  model  a 
manpower  system  as  a  markov  chain  the  following  steps  must  be 
followed. 

(1)  Divide  the  planning  horizon  to  small  intervals.  Small 
intervals  may  lead  to  more  accurate  estimates,  but  the  TPM 
of  a  very  short  time  interval  can  prevent  the  stationarity 
of  the  TPM.  Zanakis  and  Maret  [1980]  suggest  that  the  stage 
interval  should  be  determined  based  on  objectives  of  the 
study  and  the  planning  horizon.  For  long-range  plans  usually 
a  yearly  time  interval  is  selected. 


(2)  Define  states  of  the  system.  The  states  of  a  Markov 
chain  must  be  mutually  exclusive  and  exhaustive,  ie. ,  states 
should  include  each  possible  case  once  and  only  once.  For 
example,  states  of  a  Markov  chain  can  be  defined  as  levels 
of  employment  within  an  organization.  These  states  are 
exhaustive  and  mutually  exclusive  since  each  employee  has  a 
rank  and  only  one  rank. 
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(3)  Collect  data  on  the  number  of  transitions  from  each 
state  to  other  states  to  estimate  the  TPM.  In  the  absence  of 
the  transition  data  least  square  estimates  of  TPM  may  be 
obtained  using  the  historical  data  on  the  number  of 
personnel  at  each  state  at  different  periods. 

(4)  Estimate  the  TPM  and  validate  the  model.  The  sample  size 
determines  the  accuracy  of  the  TPM  estimates,  tests  are 
used  to  validate  the  model. 

A  stochastic  model  is  presented  by  Vassiliou  [1976]  for 
wastage  analysis  in  hierarchically  structured  manpower  systems. 

In  his  paper  Vassiliou  classifies  the  reasons  for  leaving  an 
organization  as  (l)  retirement  or  death,  (2)  discharge,  and  (3) 
resignation.  He  assumes  a  constant  transition  probability  for 
discharge,  retirement,  or  death  over  the  time  period.  He  then 
tests  his  assumption  and  finds  no  evidence  to  reject  it  for  the 
cases  he  studies.  The  expected  number  of  voluntary  leaves  he  is 
equal  to  the  expected  number  of  "normal*  departures,  increased  by 
the  number  of  "frustrated"  potential  departures  proportional  to 
the  cumulative  rate  of  contraction  of  the  organization  and 
decreased  by  the  number  of  "not  frustrated"  potential  departures 
proportional  to  the  cumulative  rate  of  expansion  of  the 
organization.  He  develops  necessary  equations  for  this  model  and 
examines  the  model  using  the  data  from  two  different  firms.  Then 
he  compares  his  results  with  results  obtained  from  a 
deterministic  model  and  concludes  that  his  model  describes  the 
wastage  flow  properly  and  produces  better  results. 

Young  and  Abodunde  [1979]  investigate  the  consequences  of  a 
controlled  hiring  policy  over  a  long  period  of  time  for  a 
discrete  time  model  assuming  that  the  demand  level  reaches  a 
steady  state  after  a  sufficiently  long  time.  By  assigning  costs 
to  over-  and  under-  production  and  the  use  of  mathematical 
programming  techniques,  they  develop  a  linear  programming  model 
to  obtain  an  optimal  recruitment  policy. 

Assuming  a  constant  promotion  and  wastage  level,  Abodunde 
and  McClean  [1980]  study  the  effects  of  recruitment  that  follows 
a  Poisson  process.  They  provide  expressions  to  obtain  the  number 
of  personnel  at  each  grade  level  at  each  time  period  given  the 
desired  steady  state  size  of  a  telephone  system. 

Raghavendra  [1991]  uses  a  Markov  chain  model  to  obtain 
promotion  policies  under  certain  assumptions.  He  then  translates 
promotion  policies  into  seniority  levels  and  performance  rates 
required  for  promotion. 
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Markov  Chain  Models 


The  use  of  a  Markov  chain  in  modeling  military,  government, 
and  public  agencies  manpower  supply  goes  back  as  far  as  late 
1950.  In  this  section  we  present  a  Markov  chain  model  presented 
by  Raghavendra  [1991].  This  model  is  developed  under  the 
assumption  that  no  double  promotion  or  demotion  is  permitted  and 
that  at  the  beginning  of  each  period  a  constant  proportion  of 
staff  are  to  be  promoted  to  each  level. 

Let  t  represent  the  planning  periods,  T  the  planning 
horizon,  i  and  j  the  states,  and  K  the  number  of  states  in  the 
system,  Nj^(t)  the  number  of  staff  in  grade  j  at  the  beginning  of 
the  period  t,  P^j(t)  the  probability  that  a  member  of  staff  in 
grade  i  at  the  beginning  of  period  t  will  be  in  grade  j  at  the 
beginning  of  the  next  period,  Rj(t)  the  number  of  new  recruits 
to  grade  j  during  period  t,  Wj(t)  the  wastage  factor  representing 
the  proportion  of  members  of  staff  in  grade  j  leaving  the  system 
during  period  t  due  to  retirement,  death,  resignation,  etc.,  and 
the  ratio  of  staff  promoted  to  state  j  to  the  total  staff  who 
have  joined  state  j  (promoted  and  recruited) .  Then  under  the 
above  assumptions  the  following  equations  hold: 

K 

Nj(t  +  1)='£  ;  V  j=l,2,  .  .  .  ,if.  (1) 


K 

Pij{t) +Wj^(t)=lj  V  i  =  l,2,  .  .  .  ,/C.  (2) 


Pij(t)=07  M  j  >  i+1  (3) 

Given  the  current  personnel  structure  N(l),  desired  future 
structures  N(t),  and  wastage  factors  W (t) ,  for  all  periods  under 
study,  and  given  the  proportion  of  staff  to  be  promoted  to  each 
state,  Bj,  the  promotion  and  recruitment  policies  for  K,  then  the 
higher  level  of  employment  possible  can  be  developed  using  the 
following  equations. 


From  (2)  and  (3) 


(5) 


therefore 
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^(K-1)k(1)Nx-i(1)  =Wx(2)  -Nx(l)  (l-Wx(l)  ) 

=l»x(2) 


(6) 


Since  the  proportion  of  staff  promoted  to  state  K  is  and  the 
proportion  of  staff  to  be  recruited  to  grade  K  is  1-e^  it  follows 
that 

^(/c-i)x(i)^x-i ”®x^x(2) ,  (7) 

and 

Rf^(l)=(l~eK)^K(2)  .  (8) 


From  ( 7 ) 


and 


^(X-l) (X-l)  (^)  ■ 


®X^X ( 2 ) 
W(X-1)(1) 


(9) 


^(X-1)  (X-l)  ^i)  “i"^(X-l)  ~^(X-l)x(i)  •  (^9) 

Once  estimates  of  transition  probabilities  and  number  of  people 
to  be  hired  to  state  K  at  the  beginning  of  time  period  2  are 
obtained,  similar  information  for  other  time  periods  can  easily 
be  obtained.  For  a  numerical  example  see  Raghavendra  [1991]. 

Model  Validation 


The  flow  of  personnel  is  modeled  as  a  Markov  chain  by  first, 
selecting  the  stage  interval,  second,  defining  the  states,  third, 
collecting  data,  and  forth,  estimating  the  TPM.  Once  the  model 
is  developed  it  must  be  validated.  The  TPM  is  estimated  by 
converting  the  number  of  transitions  to  row  proportions.  This 
calculation  yields  a  maximum  likelihood  estimate  of  the  true  TPM 
if  the  process  is  stationary,  ie.,  constant  over  time.  Zanakis 
and  Maret  [1980]  present  the  following  X2  tests  to  investigate 
the  stationarity  of  the  TPM  and  its  elements.  Let  s  be  the  number 
of  nonabsorbing  states,  m  be  the  total  number  of  absorbing  and 
nonabsorbing  states,  T  the  be  number  of  stages  observed,  and 
nj^  (t)  be  the  number  of  transitions  from  state  i  to  state  j 
during  period  t.  Then  nj^(t)  ,  the  total  number  of  people  available 
in  state  i  at  the  beginning  of  period  t,  p^jj(t)  ,  the  ratio  of  ^ 
people  promoted  to  state  j  from  state  i  in  time  period  t,  and  p^^- 
the  hypothesized  stationary  probability  of  an  (i,j)  transition 
can  be  described  as  follows; 
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(11) 


ID 


«i(t)  =52 

jTl 


Pij(t)-nij(t)/ni(t) 


(12) 


Pij  = 


E  "ij(t)/E  E 

t=i  t=ij=i 


(13) 


The  following  tests  may  be  used  to  test  the  stationarity  of 
the  TPM  and  its  elements  at  the  a  level  of  significance. 

(1)  The  (i,j)  transition  probability  is  constant  over  time  if 

T 

52  ni(t)  <  Xa^[T-l].  (14) 

t  =  l 


(2)  Transitions  to  a  column  state  j  are  stationary  if 
s  T 

E  'Z  <  Xa^[s{T-l)].  (15) 

i=l  t=l 


(3)  Transitions  from  a  row  state  are  stationary  if 

m  T 

E  E  [Pij-Pljit)  ]^/Pij  <  Xa^[  (in-l)  (T-l)  ]  .  (16) 

j=l  t=l 

(4)  The  entire  TPM  is  constant  over  time  if 

s  m  T 

E  E  E  <  Xa^ts(7n-l)  (T-l)  ].  (17) 

i=l  j=l  t=l 


Markov  models  are  used  to  obtain  policies  to  achieve  a 
specific  goal  (structure  of  personnel)  with  no  consideration  for 
the  cost,  possible  constraints,  or  conflicting  objectives.  In  the 
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following  section  the  goal  progranuning  approach  along  with  three 
models  developed  for  US  Army  and  the  US  Navy  are  presented. 


Linear  Goal  Programming 

Manpower  planning  can  be  formulated  as  a  linear 
goal-programming  model,  where  each  goal  to  be  attained  is 
represented  by  one  constraint.  The  model  objective  is  to  obtain 
a  compromising  solution  that  will  satisfy  all  the  goals  as  close 
as  possible,  when  the  exact  attainment  of  all  the  goals 
simultaneously  is  impossible. 

Let  Xj  denote  the  value  of  the  decision  variable.  Let 
denote  the  goal.  Let  I  denote  the  number  of  goals  to  be 
achieved,  and  J  denote  the  number  of  decision  variables.  Let 
denote  the  per  unit  contribution  of  the  decision  variable 
into  goal  i.  Let  gj  denote  the  goal  under-achievement  variables, 
and  g*  denote  the  goal  over-achievement  variables;  and  g^  are 
called  deviation  variables.  Let  wj  and  represent  the  priority 
weights  of  under-  and  over-achievement  of  goal  i.  Therefore,  a 
general  formulation  of  the  goal -programming  model  can  be 
represented  as  follows: 

I  I 

Minimize  52  53  **^1 

i=l 


j 

s.t.  J2  ^ij  Xj  +  gl  -  gl  =  Gj_;  i=l,...,m,  (19) 

j=i 


Xj,9i.gl^o. 


(20) 


The  determination  of  these  priority  weights  is  presented  in  the 
next  section. 

Many  manpower  planning  problems  can  be  formulated  as  large 
goal  programs  with  hundreds  of  constraints  representing  the 
targets  and  limitations.  Several  practical  models  have  been 
developed  for  these  problems,  such  as  the  ASCAR  and  the 
ELIM-COMPLIP  models.  These  models  and  a  Navy  model  are  presented 
in  the  following  sections. 

A  US  Navy  Model 

The  manpower  requirement  for  the  US  Navy  is  embodied  in  the 
present  and  future  needs  for  a  wide  mix  of  specialties,  with 
different  years  of  commissioned  service.  The  US  Navy  has  several 
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commissioning  sources  that  have  different  capacities  and  costs. 
Some  sources,  such  as  the  US  Navy  Academy  (USNA) ,  provide 
officers  with  wide  specialty  areas,  while  others  produce  officers 
with  a  single  specialty.  Therefore,  to  meet  the  future  needs  of 
experienced  officers  with  different  specialties,  a  proper  mix  of 
the  commissioning  sources  must  be  established.  A  flexible 
approach  is  needed  to  provide  a  fast  response  to  major  changes  in 
manpower  policies,  and  to  study  different  scenarios  of  source  mix 
and  inventory  costs. 

The  Approach.  In  this  model,  officer  inventory  within  each 
community  is  specified  by  the  number  needed  in  different  states 
at  each  time  period.  An  officer  state  is  determined  by  three 
factors:  (l)  warfare  community,  (2)  commissioning  program,  and 
(3)  years  of  commissioned  service  (YCS) .  Since  most  communities 
follow  a  common  promotion  path  in  promoting  their  officers  at 
known  YCS  experience  levels,  YCS  is  then  used  as  a  substitute  for 
grades . 

The  expected  flows  between  states  in  successive  time  periods 
is  projected,  in  a  markovian  trend,  using  transition  rates 
estimated  from  historical  data.  New  accessions  are  also  added  to 
these  flows. 

Formulation.  The  model  developed  by  Bres,  Burns,  Charnes, 
and  Cooper  [1980]  utilizes  discrete  time  periods.  The  subscripts 
used  in  this  model  are  defined  as  follows: 

i  =  Warfare  community, 
j  =  Commissioning  source, 
t  =  Time  period, 

k  =  Number  of  years  of  commissioned  service,  and 
m  =  level  of  experience  which  is  defined  by  lower  and  upper 
limits  on  YCS. 

The  constants  used  in  this  model  are  defined  as  follows: 

I  =  Number  of  communities, 

J  =  Number  of  commissioning  sources, 

K  =  Maximum  length  of  service, 

M  =  Number  of  experience  levels, 

T  =  The  planning  time  horizon, 

=  Initial  officer  inventory  in  community  i,  form  source  j, 
with  k  YCS, 

s\j(k)  =  Officer  survival  rate  for  t  time  periods  in  community  i, 
from  source  j,  with  k  YCS, 

=  Officers  strength  goal  in  community  i,  for  experience 
level  m,  in  time  period  t, 

U(t)  -  Upper  limit  on  total  officer  inventory  in  period  t, 

B(t)  =  Budget  limit  for  pay  and  allowances  in  period  t, 

=  Cost  for  pay  and  allowances  for  an  officer  in  community 
i,''from  source  j,  with  k  YCS,  in  time  period  t. 
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=  Lower  limit  for  YCS  in  experience  level  m  for  conununity  i, 
Uj^  =  Upper  limit  for  YCS  in  experience  level  m  for  conununity  i, 
XftJ  =  Upper  limit  on  the  total  number  of  officers  that  may  be 
commissioned  in  time  period  t, 

Pj(t)  *  Maximum  number  of  officers  that  can  be  commissioned  from 
source  j  in  time  period  t, 

Qj(t)  *=  Minimum  number  of  officers  that  can  be  commissioned  from 
source  j  in  time  period  t, 

Rj^(t)  »  Maximum  number  of  newly  commissioned  officers  that  can  be 
trained  in  community  i,  in  time  period  t, 

P^j(t)  =  Maximum  allowable  number  of  officers  commissioned  from 
source  j  to  community  i,  for  time  period  t, 

Q^j(t)  =  Minimum  allowable  number  of  officers  commissioned  from 
source  j  to  community  i,  for  time  period  t, 

=  Weight  given  to  over-achievement  of  officer  strength 
goal  for  community  i,  experience  level  m,  time  period  t, 

W]j„(t)  =  Weight  given  to  under-achievement  of  officer  strength 
goal  for  community  i,  experience  level  m,  time  period  t,  and 
w~(t)  =  Weight  given  to  negative  deviation  from  total  officer 
strength  limit,  time  period  t. 

The  decision  variables  used  in  this  model  are  defined  as  follows: 

yijk(^)  -  Officers  inventory  in  community  i,  from  source  j,  with 
k  YCS,  at  the  beginning  of  time  period  t. 

Xj^^(t)  =  Accessions  to  community  i  from  source  j  in  time  period 

t. 

g'±m(t)  =  Goal  under-achievement  for  community  i,  experience  level 
m.  in  time  period  t. 

g±„(t)  =  Goal  over-achievement  for  community  i,  experience  level 
m,  in  time  period  t. 

g~(t)  =  Negative  deviation  from  the  total  officer  inventory  upper 
limit  in  period  t. 

The  constraints  that  define  the  relation  between  officer 
inventories  and  goals,  for  each  community  at  different  experience 
levels,  are  expressed  as: 

J  “im 

E  T.  9lm(^)=  Gjjnit)  ,  (21) 

where  inventories  can  be  expressed  in  terms  of  beginning 
inventories  as: 

yij;c(t)  =  for  2<t<k,  (22) 


and  in  terms  of  subsequent  accessions  as: 

The  limitation  on  total  officer  inventories  is  expressed  by: 
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for  t>k . 


(23) 


yijk(t)=  s}j(0)xj_j(t-k) 

I  J  K 

53  V  yij7c(t)+  srla(t)=  U(t)  forall  t.  (24) 

k=l 


The  budget  limits  are  expressed  by: 

I  J  K 

E  E  E  ^ljk(t)yijkit)  <  Bit)  V  t.  (25) 

i=l  j=l  *=1 

Restrictions  on  the  officer  accessions  from  various  warfare 
communities  are  addressed  in  several  constraints.  The  limitation 
on  the  total  number  of  the  officer  accessions  in  each  period  is 
expressed  by: 

Xijit)  <  Pit)  V  t.  (26) 

i=i  j=i 

In  each  community,  the  training  capacity  limits  for  newly 
commissioned  officers  are  expressed  by: 

J 

Yxijit)  <  Riit)  Vi,t.  (27) 

J=1 

The  operating  upper  and  lower  limits  for  the  commissioning 
sources  are  expressed  by: 

I 

Qjit)  <  53xij(t)  <Pjit)  Vi,t.  (28) 

i=l 

The  upper  and  lower  limits  on  the  distribution  of  newly 
commissioned  officers  from  each  source  to  each  community  are 
expressed  by: 

Oij(t)  <  Xij(t)  <  Pij(t)  V  t.  (29) 

Let  ot(t)  ,  P(t),  '1  (t)  ,  and  S  (t)  denote  the  minimal  and  maximal 
proportions  of  change  allowed  between  adjacent  time  periods.  The 
proportional  changes  allowed  for  inputs  to  each  community  are 
expressed  as  follows: 

for  the  accessions  to  community  i.  The  proportional  changes 
allowed  for  outputs  from  each  source  are  expressed  as  follows: 
for  the  output  of  source  j.  Finally,  a  nonnegativity  condition  is 
imposed . 
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J  J  J 

“(t)  i  ^(t)Vxij(t)  V  t  and  (30) 

j=i  jTi  fTi 

II  I 

T(t)  53  XijCt)  <  5^  5  tf(t)  ^ij(t)  t  and  (31) 

i=l  i  =  l 

The  model  objective  is  to  find  an  officer  accession  and  an 
arrangement  plan  that  minimizes  weighted  deviations  from  officer 
strength  goals.  This  can  be  expressed  as  follows: 


TIM  T 

Minimize  E  (‘^in(t)gr^(t) +>^ia(t)gria(t) )  +J^v-(t)g~(t), 

t=li=lm=l  t=l  (32) 

subject  to  constraints  (21) -(31). 

As  an  illustration,  the  model  was  applied  to  four 
communities  (Unrestricted  Line  officers)  for  the  first  ten  years 
of  commissioned  service  omitting  the  budget  constraint.  It 
provided  a  satisfactory  accession  plan  with  minimal  deviations 
for  the  high  priority  requirements.  For  more  details  see  Bres, 
Burns,  Charnes  and  Cooper  [1980].  The  model  can  also  be  used  to 
provide  the  impact  on  budgets  if  the  appropriate  results  are 
substituted  into  constraint  (25)  to  study  funds-flow  consequences 
of  these  plans. 

The  ASCAR  Model 


The  Accession  Supply  Costing  and  Requirement  (ASCAR)  model 
is  a  goal  programming  based  model.  It  was  initially  developed  by 
General  Electric  for  the  Congressional  Budget  office  in  the  late 
seventies.  The  model  was  then  revised  and  used  by  the  office  of 
the  Assistant  Secretary  of  Defense.  The  main  objective  of  the 
ASCAR  project  was  to  provide  the  capability  to  investigate  the 
effects  of  changing  manpower  policies  on  personnel  requirements, 
Tiilitary  qualifications,  and  the  associated  costs. 

The  ASCAR  model  consists  of  data  management  programs,  a  data 
base,  personnel  flow  simulation  routines,  a  cost  routine,  and  a 
report  routine.  It  follows  a  five-step  process  to  evaluate  the 
annual  accessions  necessary  to  meet  the  strength  and  quality 
requirements,  and  to  calculate  the  associated  costs.  These  steps 
are  summarized  in  the  following: 

•  Initially,  the  historical  data  are  analyzed  to  develop 
flow  rates  and  starting  personnel  levels  for  the  simulated 
period. 
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•  These  rates  and  levels  are  then  used  in  the  personnel  flow 
simulation  routines  to  compute  losses,  to  forecast  the 
supply  of  new  recruits  for  each  simulated  year,  and  to 
simulate  the  flow  of  new  recruits  in  order  to  calculate  the 
required  accessions. 

•  The  supply  of  new  recruits  is  categorized  by  sixty 
qualitative  factors  such  as  educational  level  and  mental 
category . 

•  The  goal  programming  procedure  optimizes  the  different 
category  mix  of  new  recruits  to  match  the  required 
man-strength  and  characteristics  while  satisfying  the 
specified  constraints. 

•  The  cost  module  produces  the  cost  estimates  for  the 
analyzed  policy  using  the  cost  factors  and  the  projected 
strength  determined  in  the  previous  steps. 

The  process  is  then  repeated  for  each  additional  year.  More 
details  are  provided  in  Collins,  Gass,  and  Rosendahl  [1983]. 

The  ELIM-COMPLIP  Model 

The  Enlisted  Loss  Inventory  Model-Computation  of  Manpower 
Programs  Using  Linear  Programming  (ELIM-COMPLIP)  is  a  goal 
programming  based  forecasting  system  developed  by  General 
Research  Corporation  for  the  Department  of  the  Army.  It  is  used 
for  manpower  planning,  budgeting,  and  personnel  policy 
formulation.  The  system  is  described  in  details  by  Holz  [1980]. 

In  the  late  sixties  the  Assistant  Secretary  of  the  Army, 
William  K.  Brehm,  was  addressing  several  ''what  if"  questions 
related  to  different  manpower  policies  for  Vietnam.  Each 
alternative  required  approximately  8-20  computation  hours. 
Therefore,  the  General  Research  Corporation  was  asked  to  automate 
these  computations.  One  of  the  two  systems  developed  was  the 
COMPLIP  optimization  system.  The  main  objective  of  COMPLIP  is  to 
minimize  the  weighted  sum  of  the  deviations  of  the  actual 
man-strength  from  the  required  goals.  It  provided  the  fast 
responsiveness  needed  to  generate  accurate  manning  plans  for 
models  with  conflicting  constraints  that  were  difficult  to 
approximate  by  manual  calculations.  The  second  was  a  simulation 
of  the  manual  system.  It  as  never  used  after  the  demonstration 
of  the  COMPLIP  system.  COMPLIP  provided  upper  and  lower  bounds  on 
draft  calls  that  were  very  useful  in  preventing  fluctuations  to 
this  sensitive  political  issue. 

In  the  fiscal  year  1972,  the  Army  was  phasing  down  from 
Vietnam,  when  the  Congress  passed  an  authorization  bill  to  reduce 
the  Army's  strength  by  50,000  for  that  yea',  which  significantly 
accelerated  the  phase-down  plan.  Several  drastic  policies  were 
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adopted  that  disrupted  personnel  management  and  caused  the 
strength  to  fall  below  authorized  level.  The  Army  ended  the  year 
50,000  men  below  the  desired  strength.  As  a  response  to  this 
crisis,  the  General  Research  Corporation  developed  the  Enlisted 
Loss  Inventory  Model  (ELIM)  using  a  new  loss  projection  method. 
The  new  method  starts  with  the  most  recent  information  and 
preserves  the  relations  between  end  strength,  gains,  and  losses 
to  avoid  the  risks  created  with  the  previously  used  process. 

The  ELIM-COMPLIP  system  is  self -correcting  with  respect  to 
projection  errors,  and  provides  accurate  strength  forecasts.  It 
includes  five  modules:  data  processor,  rate  generator,  inventory 
projection,  optimization,  and  a  report  generator  module.  The 
primary  task  of  the  rate  generator  is  to  exponentially  smooth  the 
time  series  of  historical  data  to  project  loss  rates.  These  rates 
are  used  in  the  inventory  projection  module  to  compute  retention 
rates  that  develop  the  goal  programming  coefficients.  For  more 
details  see  the  description  and  the  system  schematic  in  Holz 
[1980] . 


Linear  goal-programming  has  proved  to  be  very  effective  in 
solving  multi-criteria  problems,  such  as  manpower  planning 
problems.  Usually,  the  objective  function  (18)  of  a  general 
linear  goal-programming  model  (18) -(20)  contains  several 
thousands  of  deviation  variables.  Most  problems  have  specific 
properties  that  necessitate  proper  determination  of  weights,  wj 
and  w^,  associated  with  the  deviation  variables  to  reflect  the 
correct  goal  priority.  In  general,  direct  weight  determination 
for  small  problems  with  related  goals  is  not  a  hard  task. 

However,  direct  weight  determination  for  thousands  of  dissimilar 
goals  is  not  viable,  as  it  resembles  comparing  ''oranges  to 
apples".  In  this  section  we  represent  a  weight  determination 
approach  proposed  by  Gass  [1986]  using  Saaty's  analytical 
hierarchy  process  (AHP) . 

Consider  the  grade-skill  goal  constraints  (21) ,  budget 
constraints  (25) ,  and  promotion  goal  constraints  (30) .  Meeting 
any  of  these  constraints  has  a  completely  different  meaning  than 
meeting  the  other  two.  Also,  meeting  a  grade-skill  goal  in  one 
year  may  have  a  much  different  priority  than  the  same  goal  in  a 
different  year.  Therefore,  planners  must  be  able  to  designate  the 
objective  goal  weights  that  cause  a  proper  transition  to  be  made 
over  the  planning  horizon  to  reach  specified  requirements  at 
certain  years.  Undoubtedly,  the  exact  matching  of  man-strength 
requirements  that  fluctuate  over  the  intended  period  is  not 
achievable.  However,  a  compromising  solution  can  be  accomplished 
where  the  goal  weights  will  determine  the  years  to  be  closely 
matched,  namely,  the  solution  profile. 


Gass  [1986]  developed  an  automated  routine  for  calculating 
thousands  of  weights  necessary  to  determine  an  acceptable 
solution  for  the  linear  goal  program  (18) -(20). 

A  grade-skill  model  with  a  seven  year  planning  horizon  is 
•sed  to  illustrate  Gass'  method.  In  this  model  the  goals  are 
defined  by  (1)  total  personnel  at  the  end  of  each  year,  (2) 
grade-skill  qualifications,  (3)  promotion  goals  by  grade-skill, 
and  (4)  loss  goals  by  grade-skill.  These  goals  can  be  represented 
by  a  hierarchical  structure  of  factors  illustrated  in  Figured). 
Level  1  represents  the  model  main  objective.  Level  2  contains  the 
major  requirements  (year  end-strength)  that  influence  the  model. 
Level  3  contains  the  grade-skill,  promotion,  loss,  and  gain 
targets  that  influence  level  2. 


1986) 


The  Gass  method  establishes  the  overall  priorities  in  the 
main  objective  from  the  hierarchical  factor  priorities  by 
systematically  comparing  the  factors  within  each  level.  In  each 
level,  the  weights  that  reflect  the  factor  priorities  are 
obtained  from  a  comparison  matrix  associated  with  the  level.  The 
matrix  elements  have  values  1  to  9  and  their  reciprocals,  and 
represent  the  importance  of  meeting  each  factor.  A  typical 
comparison  matrix  corresponding  to  level  2  of  the  above  example 
is : 
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1  5  3  1  3  3  3 

1/5  1  1/5  1/3  1/5  1/7  1/9 

1/3  5  1  7  7  9  9 

131/7  1  3  5  7 

1/3  7  1/9  1/5  1/3  1  5 

1/3  9  1/9  1/7  1/3  1/5  1 

where  Aj  ^=1  means  that  year  1  and  year  4  are  equally  important, 
and  Ay  2=^  means  that  year  7  is  significantly  more  important  than 
year  2',  with  respect  to  the  model  main  objective.  The  Gass  method 
determines  the  solution  to  the  system  of  equations  Aw=  Aw,  which 
is  the  well-known  eigenvalue  problem.  Following  the  AHP  theory, 
for  A  equal  to  the  largest  eigenvalue,  the  normalized  w^ 
components  are  then  interpreted  as  weights  that  represent  the 
importance  of  each  factor.  For  the  above  matrix  we  have  the 
following  weights  for  level  2: 

Year  1  2  3  4  5  6  7 

Weights  0.25  0.03  0.38  0.18  0.07  0.08  0.01 

These  weights  indicate  that  the  accomplishment  of  year  3's  goal 
has  the  highest  priority  (0.38),  then  year  I’s  goal  and  so  on. 

A  similar  analysis  is  carried  out  to  determine  the  level  3 
factor  priorities  for  each  year.  The  composite  goal  priorities 
for  the  level  3  items  are  developed  by  multiplying  the  year's 
priorities  by  the  corresponding  factor  priority,  and  then  adding. 
A  detailed  example  is  described  in  Gass  [1986]. 

Finally,  the  AHP  priorities  are  converted  to  goal  weights 
for  each  target-year  combination.  This  is  accomplished  using 
closeness  factors.  Each  target-year  combination  is  characterized 
by  an  ordered  couple  c(i,j) ,  in  this  example  i=l,...,4  (targets) 
and  j=l,...,7  (years),  which  adjusts  the  composite  weight  for 
target  i  to  account  for  the  priority  of  year  j.  Thus,  weight  for 
any  goal  i  and  year  j  is  determined  by 
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(33) 


w^j=Si-c{i,  j)  (  ) 


where  S*  is  the  top  of  the  scale  for  goal  i  and  sf  is  the  bottom 
of  the  scale  for  target  i.  A  broad  scale  (e.g.  0  to  1000)  is  used 
so  that  the  optimization  algorithm  would  be  able  to  differentiate 
properly  between  the  goals. 

The  AHP  theory  can  be  extended  to  determine  the  skill 
weights  for  each  grade,  in  order  to  set  priorities  between  the 
skills  within  a  grade  for  each  year.  An  example  of  weights 
indexed  by  grade  and  skill  is  given  in  Gass  [1986]. 

Markovian  Goal  Programming 

Even  though  Markov  models  are  successfully  used  to  estimate 
the  future  manpower  plan  or  to  develop  policies  to  obtain  a 
desired  future  personnel  structure,  they  cannot  consider  cost, 
constraints,  or  conflicting  objectives.  Zanakis  and  Maret  [1981] 
use  a  combination  of  the  Markovian  and  LGP  techniques  to  obtain  a 
manpower  plan.  As  previously  stated  they  use  historical  data  to 
estimate  TPM  and  desired  structure  of  the  personnel  system.  Once 
the  number  of  personnel  is  determined  they  develop  a  goal 
programming  model  to  satisfy  the  organizational  goals  such  as: 
keep  the  cost  at  the  lowest  possible  level,  hire  a  specific 
number  of  recent  college  graduates,  keep  the  ratio  of  number  of 
personnel  at  a  specific  level  of  employment  to  other  levels 
constant,  and  other  similar  goals.  For  details  and  a  numerical 
example  see  Zanakis  and  Maret  [1981]. 

Assuming  that  costs  can  be  assigned  to  under  and  over 
production  and  given  the  time  demand  level  will  reach  steady 
state.  Young  and  Abodunde  [1979]  present  a  formulation  for  a 
discrete-time  Markovian  goal  programming  model. 

Let  G(t)  denote  the  target  demand  at  time  period  t,  and 
C~(t)  and  C*(t)  denote  the  costs  ascribed  to  the  goal  under-  and 
over-achievement.  Let  R(t)  denote  the  total  recruit  inventory  at 
time  period  t,  and  the  number  of  units  surviving  the  first  s 
years  of  service.  Let  P  denote  the  transition  probability  matrix, 
Hq  denote  the  recruitment  vector,  and  (l~a)  denote  the  failure 
rate.  Then,  the  linear  goal  programming  model  (18) -(20)  can  be 
stated  as  follows: 

T 

Minimize  ^  {  C*  it)  g*  [t)  +  C  {t)  g~  it)  )  (34) 

C»1 
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t 

S.t.  ^3  Rt*l-a*  9~(t)-  g*(t)  =  G{t)  ;  t  =  l,...,T, 

3  =  1 


(35) 


(36) 


Notice  the  resemblance  of  (35)  and  (36)  to  (21)  and  (22) . 

Other  constraints  such  as  limits  on  the  number  of  recruits 
for  some  or  all  of  the  periods  under  study  can  be  added.  The 
technique  is  applied  to  a  case  study  based  on  the  Irish  Telephone 
System.  Numerical  results  are  presented  in  Young  and  Abodunde 
[1979]  . 

Conclusions 

Markov  models  are  valuable  tools  in  manpower  planning. 
Stochastic  and  deterministic  models  are  developed  and 
successfully  used  to  predict  future  manpower  structures,  to 
derive  policies  to  achieve  a  specific  manpower  structure  in  the 
near  and  far  future,  and  to  provide  insight  to  *what  if'» 
questions.  Historical  data  are  used  to  estimate  the  TPM.  The 
estimated  TPM  is  not  very  accurate,  so  the  results  from  Markov 
models  must  be  used  cautiously.  Deterministic  models  are 
preferred  by  some  practitioners  since  they  believe  that  forecasts 
are  not  accurate  enough  to  justify  the  complexity  of  the 
stochastic  models.  Even  though,  Markov  models  are  very  helpful  in 
manpower  planning  they  do  not  take  into  account  costs  or 
conflicting  objectives. 

Manpower  planning  is  closely  linked  to  linear  programming 
theory.  Manpower  planning  problems  can  be  represented  as  a  linear 
goal-programs,  where  limitations  and  conflicting  goals  are 
formulated  as  constraints  with  objectives  to  minimize  the 
deviations  from  these  limitations.  Several  practical  models  have 
used  this  approach  for  solving  manpower  problems.  This  approach 
has  proven  to  be  very  efficient  in  providing  satisfactory 
results.  In  goal  programming  model,  goal  priorities  have  a 
significant  role  in  determining  how  closely  each  goal  is  matched, 
and  consequently,  the  solution  profile.  One  of  the  suggested 
techniques  for  determining  the  objective  weights,  that  reflect 
the  priorities,  is  to  systematically  compare  the  main  factors  at 
different  levels  of  the  problem.  Most  of  the  data  used  in  these 
goal  programs  rest  on  assumptions  and  predictions  about  several 
human  and  economic  factors,  e.g.  officer  survival  rate.  To 
properly  incorporate  these  nondeterministic  factors,  Markov 
models  are  used  to  quantify  these  predictions  to  be  used  as 
goal-program  constants. 
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The  Prototype  Model 


The  model  described  below  is  designed  so  that  two  reports 
can  be  generated  and  used  for  decision  making.  The  inputs  to  the 
model  are  the  following; 

(i)  the  force  structure, 

(ii)  inventory  targets, 

(iii)  promotion  targets,  and 

(iv)  transition  rules  for  personnel  who  move  through  the  system. 
The  output  tables  take  the  following  form: 


Personnel  Inventory  Report 


Time  Grade  Policy  Target 
Group  (Input) 


Inventory  Deviation 
Produced 
By  The 
Model 


Deviation 


Personnel  Promotion  Report 


Target  %  promoted 
at  this 
(Time, Grade) 
combination 
(input) 


%  Promotions 
produced  by 
the 
model 
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Graphical  Display 


policy  group  i 


Subscripts 

t-  denotes  the  time  period  (years) 
g-  denotes  the  grade 

r-  denotes  the  location  (r=l  conUS,  r=2  not  conUS) 
i-  denotes  a  policy  group 

m-  denotes  the  number  of  years  in  policy  group  i 
It  is  assumed  that  all  personnel  can  be  characterized  by  a  5- 
tuple  given  by  the  5  subscripts.  We  will  refer  to  a  5-tuple 
(t,g,r,i,m)  as  a  state  from  which  a  transition  can  be  made  to 
another  state. 
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Graphically  this  may  be  viewed  as  follows: 


The  rules  for  the  possible  transitions  are  input  to  the  model. 
Node  Names 


(t,g,r,i,m)  -  state 

ACC  -  denotes  the  origin  node  for  all  accessions 

SEP  -  denotes  a  destination  node  for  all  separations 

BAL  -  denotes  a  balance  node  to  balance  the  supply  and 

demand  in  the  model.  This  is  required  for  a  true  network  model  of 
the  form  AX=r  where  A  is  a  node-arc  incidence  matrix. 


Arc  Types 


(t , g, r, i, m; t+1, g, f , i, m)  -  arcs  which  corresponding  to  moving  from 
state  (t,g,r,i,m)  to  state  (t+1,  g,  f ,  i, m) 

(ACC;t, 1, r, 1, 1)  arc  from  the  node  ACC  to  state  (t,l,r,l,l). 

(t , g, r, i, m; SEP)  arc  from  the  node  (t,g,r,i,m)  to  the  separation 
node . 

(H,  g, r, i, m;BAL)  arc  from  state  (H,g,r,i,m)  to  the  node  BAL. 

(SEP; BAL)  arc  from  node  SEP  to  node  BAL. 
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Figure  3 .  General  structure 


Network  structure 


Constants  and  Sets 


H-  planning  horizon  (t=l,...,H) 

T(g,r,i,m)  -  the  set  of  nodes  (g,f,i,m)  for  which  there  exists  an 
arc  (t,g,r,i,m;t+l,g,f ,i,m)  for  t=l,...,H-l.  This  is  sometimes 
called  the  to  set  or  the  after  set. 

F(g,r,i,m)  -  the  set  of  nodes  (t-l,g,r, I,m)  for  which  there 
exists  an  arc  (t~l,g,f ,1 ,m;t,g,r, i,m)  for  t=2,...,H.  This  is 
sometimes  called  the  from  set  or  the  before  set. 

RHS(g,r,i,m)  -  niimber  of  people  in  the  state  (l,g,r,i,m) 

ES(t,g)  -  force  structure  (people) 

N(t,g,i)  -  inventory  targets  (people) 

PROM(t,g)  -  Promotion  Targets  (%) 

W2,  W3,  W4  -  weights  for  deviations  from  targets. 

Decision  Variables 

X(t,g,r, i,m;t+l,g,f ,i ,m)  -  number  of  people  on  the  arc 
(t,g,r,i,m;t+l,g,f ,t,m) 

I(t,g,i)  -  number  of  people  who  leave  state  (t,g,i) 

P(t,g)  -  number  of  people  who  leave  state  (t,g)  and  are  promoted 
A(t,r,l,l)  -  number  of  people  on  arc  (ACC,t, l,r, 1, 1) 

B(g,r,i,m)  -  number  of  people  on  arc  (H,g,r, i,m;BAL) 

C  -  number  of  people  on  arc  (BAL,ACC) 

D  -  number  of  people  on  arc  (SEP,BAL) 

S(t,g,r,i,m)  -  number  of  people  on  arc  (t,g, i,m;SEP) 

IO(t,g,i)  -  number  of  people  over  the  desired  inventory  level 
IU(t,g,i)  -  number  of  people  under  the  desired  inventory  level 
PO(t,g)  -  number  of  people  over  the  desired  promotion  %  goal 
PU(t,g)  -  number  of  people  under  the  desired  promotion  %  goal 
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Constraint B 


Conaervation  of  Flow  (t .a.  i.in)  =  (  1, 1, 1. H 


*  SU.l.z.l. 

(5.r,T,S)  eT(2,r,  i,»)  . 

=  RUSd.r.l.l)  V  Ir.i.rn) 


1)  -  Ad.r.l.l) 


|ifM=l 
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Conservation  of  Flow  t=l,  (g, i,m) / (1, 1, 1) 

xn,  g,  r,  +  S{l,g,  r,  i,m) 

(•y,T',T,W)eT(g,  r,  l,m) 

=  RHS  (g,  r,  i,  m)  V  (g>l,r,i,m) 
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Conservation  of  Flow  l<t<H  ,  g=l 

X(t,  1,  r,  i,m;  r,T,7ff)  +  S(t,  1,  r,  i,/n) 

r,  i.m) 

if  (t-l,  1,  r,T,7JJ;  t,  1,  r,  i,  m) 

(g,  r,  l,m)€f  (1,  r,  i.ai) 

-A{t,r,i,m)  =  0  V  1,  r,  i, jn) 
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Conservation  of  Flow  l<t<H,  (g,  i,  m)  (1/ 1/ 1) 

X(t,g,r,  i,m;t+l,g,r,r,m 

Cy,T,1',W)€T(g,  r,  l,m) 

5^  X{t-l,g,r,T:,m:t,g,r,i,m)  =  0 

C7,T,T,lir)*f<l,  r,  l,m) 


+  S(t,  1,  r,  i,;n) 
V  {\<t<H,g,r, 
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Conservation  of  Flow  t=H,  (g, i,m) = (1, 1, 1) 

S(H,l,r,l,l) +B(l,r,l,l)  -A(H,r,l,l)  =0 


Figure  9.  Conservation  of  flow  t=H,  (g,i,in)  =  (l,l,l) 


Conservation  of  Flow  t=H,  (g, i,m) ^ (1, 1, 1) 

S(H,  g,  r,  i,m)  -  X(H-1  ,'§,Y,'x,TR;  H,  g,  r,  i,  m) 

fy,  T,  T,  jr)  €  f  ( gr,  r ,  i ,  a ) 

+  B{g,r,i,m)  =  0,  V  (g,r,i,m) 


SEP  SEP 

Figure  11.  Conservation  of  flow  at  ACC 
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Conservation  of  Flow  at  SEP 

D  -  S{t,g,r,i,m)  =0 

t,g,  r,l,a) 

Inventory  Constraints  t<H 

r,a  f7,T,T,W)eT(g,r,  1,«) 

+Y^S(t,g,r,i,m)  -Ht,g,i)  =  0,  V  (t<H,g,i) 

r,  m 


I' 


I  all  ^  Kt.gj^ 


Figure  12.  Inventory  constraints  t<H 
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Inventory  Constraints  t=H 

r,  i,m)  g,  r,  i,m) 


r,m 


-  I  (H,  g,  i)  =  0,  V  (gr,  i) 


Figure  13.  Inventory  constraints  t=H 
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Force  Structure  Constraints 


Y,I{t,g,i)  =  ES{t,g) ,  V  (t,g) 

I 

Promotion  Accounting  Constraints  t<H, g<G 

52  52  X(t,g,r,i,in;t+l,g+l,r,i:,m) 

(r,  i,/n)  l,m) 

-  PU,g)  =  0  V  lt<H,g<G) 


Figure  14.  Promotion  accounting  constraints  t<H,  g<G 
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Deviation  From  Inventory  Goal 

I{t,g,i)  +  IU(t,g,i)  -  lO(t,g,i)  =H{t,g,i),  V  (t,g,i) 

Deviation  From  Promotion  Goal 

P(t,g)  ^  PU(t,g)  -PO(t,g)  =  PROM[t ,g)ES{t ,g)  ,  V  {t<H,g) 

Objective  Function 

minimize  ^  {W^IU(t,g,i)  +  W2lO(t,g,i)) 

+  j;  (P/3PU(t,g)  +  f/4P0(t,g)) 

(t,g) 


Constraint  Layout 


Name 

Size 

States 
t,g,r, i,m 

ST 

H.G.R.I.M 

ACC 

ACC 

1 

SEP 

SEP 

1 

BAL 

BAL 

1 

Inventory 

INV 

H.G. I 

Force 

FOR 

H.G 

Promotion 

PRO 

(H-1) . (G-1) 

Inventory  Goal 

IG 

H.G. I 

Promotion  Goal 

PG 

(H-1) . (G-1) 

Total  Constraints  =  H.G.R.I.M  +  2.H.G.I  +  2 . (H-1) . (G-1) +  H.G  +  3 
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Appendix  A 

Model  Generator  and  Test  Case 
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